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The back-to-back J/ip and tt associated production at e''"e“ colliders is proposed to detect the 
gluon transverse momentum dependent(TMD) fragmentation functions. TMD factorization is as¬ 
sumed for this process. With spinless pion, unpolarized and linearly polarized gluon TMD fragmen¬ 
tation functions can be defined. It is found at parton level the hadronic tensor can be described by 
four structure functions. As a result, there are three independent angular distributions, of which a 
cos 2(j) azimuthal asymmetry is sensitive to the linearly polarized gluon fragmentation function. 


I. INTRODUCTION 

The transverse momentum dependent fragmentation function(TMDFF) is an important component of transverse 
momentum dependent(TMD) factorization [U[5]. It tells us how the hadron in a jet is affected by the transverse motion 
and polarizations of the fragmenting parton. Up to now, the quark TMDFFs have been studied very thoroughly(for a 
review, see[3]). But for gluon TMDFFs, the study is very limited. In this paper, we point out these functions can be 
extracted from quarkonium and pion associated production at e+e” colliders, i.e. e’*' -|- e“ — )■ J / ifi -\- tt -\- X . We will 
use nonrelativistic QCD(NRQCD)[3] to describe the production of J/ip- In this framework, a pair of heavy quarks is 
first produced from the hard interaction and then evolves into a quarkonium according to NRQCD. Thus at leading 
order of Ug, it must be a gluon to fragment into the final pion. If we demand J/if and tt are nearly back-to-back, 
their relative transverse momentum distribution will be very sensitive to the transverse motion of the fragmenting 
gluon, and the cross section should be described by gluon TMDFFs. Since the initial state is colorless, there is no 
interference between initial and final states and then all soft divergences can be absorbed into the fragmentation 
function and NRQCD matrix elements and a proper soft factor. In this sense we expect TMD factorization to hold 
for this process. A potential problem is whether the fragmentation function is general, or is the same as that derived 
from Semi-inclusive deep inelastic scattering(SIDIS). This can be checked by one-loop calculation. But here we just 
confine ourselves to the tree level and discuss how much information can be extracted assuming the factorization. 
The organization of this paper is as follows: In Sec.II the notations and kinematics are introduced; In Sec.Ill the 
formalism and the derivation of three independent angular distributions are given; Sec.IV includes our main result 
and a discussion; Sec.V is the summary. 


II. KINEMATICS 


The process we study is 


e+{l') + e-{l) -7 J/^{Pi) + tt{P 2 ) + A, (1) 

where V, I, Pi, P 2 are the momenta for each particle, respectively. For this process, it is convenient to work in the 
hadron frame[S]. This is a frame in the center of mass system(CMS) of leptons and with -|-z-axis along the three 
momentum of pion P 2 . The -|-a:-axis can be an arbitrary fixed direction denoted by n, which is perpendicular to P 2 . 
In this frame the momentum of lepton I is along the direction (0, tIj), as shown in Figj^a). In our interested region, 
J/ijr and tt are nearly back-to-back, i.e. < ^qcd- Then only the transverse direction of Pi is relevant. The 

azimuthal angle between l± and Pi_l is defined as </), as shown in Figj^b). Equivalently, (j) is the angle between the 
plane expanded by {P2, 1 ) and that by {P2, Pi). 

In our interested region, the invariant mass of leptons = {I + l')^ is much higher than the typical 

nonperturbative scale ^qcd- region, the mass of pion can be ignored, i.e. P 2 = 0. Since the quarkonium mass 

Mj ~ 2Mq is also a hard scale, we keep the mass of J ftp explicit in our calculation, and define 

_ Mj _ 2 Mq 

Now the TMD factorization is expected to be applicable in the region P^j^ ^ and P^^^ <C Mj. 
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FIG. 1. The scattering angles defined in hadron frame, where the momentum of pion P 2 is along + 2 -axis. I and Pi are the 
momenta of electron and J/tp, respectively. 


The differential cross section can be written as: 


da 

dz2dQdzid‘^ Pi± 


2 2 

32Tr'^Q‘^^y zf — 




(3) 


where 12 = {9, ip) is the solid angle of electron, e, eg are the electric charges of electron and heavy quark respectively, 
and zi, Z 2 are the energy fractions of J/ip and tt in CMS frame, that is, 


■zi 


Pl-q 


Z2 = 


P 2 ■ q 


(4) 


The leptonic and hadronic tensors are the standard ones, 

Lfiv — I ■ I'g^v) ~ , 

^Y.{Q\fmJli’{Pi)AP2)X){J/iP{P^)^{P2)X\3^^mQ)5\q - Pi - P 2 - Px), 


X 


(5) 


where = ip^^ip is electromagnetic current and A is the helicity of electron. 

For the calculation, it is convenient to define the transverse direction through two light-like vectors: P 2 and 
q = q — 2 ^P 2 - The transverse metric and anti-symmetric tensor are defined as: 


= g^'' — n^n" — n" n 
1 

iq 








n„n 


p' ? 


VP 2 -q' 


y/P 2 ■ q 


^2 ’> 


( 6 ) 


with £0123 _ light-cone coordinates the -I-, — components of a vector are a~^ = n ■ a, a = n ■ a, and the 

transverse component is = g^ a^. 


III. FORMALISM AND STRUCTURE FUNCTIONS 


At leading power(or twist) of Pi±IQ and Pij_/Mq, the cross section can be factorized into the product of a hard 
kernel and a fragmentation function for the pion production, as shown in Figj^ The fragmentation function is 
transverse momentum dependent and can be defined as [6]: 




(7) 


where is the gluon field strength tensor, a is the color index. Here we have assumed the hadron h is moving 
along -|-z-axis, then the large component of its momentum is P^. The parton momentum fraction is 1/z = k~^/Pj^, 
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FIG. 2. The leading region for Jip, tt back-to-back production. The central bubbles represent the hard scatterings. 


0 < z < 1. We have suppressed the gauge links, which are irrelevant to our discussions here. From the definition, G 
corresponds to the unpolarized fragmenting gluon and H corresponds to the linearly polarized gluon. Since the final 
hadron is spinless or unpolarized, there are only these two fragmentation functions at leading twist. 

For the quarkonium production, we use NRQCD|4]. The heavy quark pair QQ is first produced from the hard 
interaction, then transforms into a quarkonium. The QQ has a small relative velocity v in the rest frame of the quark 
pair. NRQCD gives a consistent power expansion according to v. To a certain power of v, only finite number of 
NRQCD matrix elements contribute. For the process considering now, the quark pair from hard interaction must be 
in color octet. To leading order of v, there are three types of NRQCD operators contributing to J/4’ production: 
O8(1^0),O8f^l), 08^^7)111151. The quantum number represents the angular momentums of the quark pair 

in the rest frame. Here we use the covariant formalism in 1^ to project the quark pair to definite partial waves. For 
details please see l9], here we just list the final results. 

According to the collinear power counting rule [HU ITT] . at leading twist the fragmenting gluon is collinear to the 
final pion, i.e. = {k'^, k~, kgj_) ~ (5(1, A^, A) with A ~ ^qcdIQ ^ 1- This results in the leading region shown 
in Fig[^ Contrary to the collinear factorization, now the transverse momentum kg± cannot be ignored in the hard 
part. Actually, this momentum is contained in the delta function for momentum conservation, i.e. 5^{Pi± + kg±). 
Since this delta function has already been of leading power, kg± can be ignored in other parts of the hard scattering 
amplitude. Following the procedure in e.g. the resulted hadronic tensor is 


W^'' = 


2z2 


(1 - t2)Q2M, 


-S{zi - zt) J + fcijTWphgf G( 2 ,fc 


2k1kl - • kj_ 


2M1 


H{z,kl)]. (8) 


where fivpT are Lorentz indices, as shown in Figj^ Notice that in this formula the energy fraction of J/'i/' 2:1 is totally 
fixed to be zj = (1 -I- r^)/2. This is a result of Pi± —>■ 0. Moreover, the momentum fraction z is determined by zi 
and Z 2 definitely. With zi = zj, z and Z 2 satisfy a very simple relation: 


1 _ 1-t2 
Z 2z2 


(9) 


The tensor is the module of the hard amplitude, 

^pr = E — A>;,{A';nosC^+^Lj)). 

LSJ 


( 10 ) 


where we have summed over all possible partial waves with S' < 1 and T < 1; is the hard partial wave amplitude, 
but without the polarization vectors for fragmentation gluons. The average over the color and polarizations of the 
quark pair is necessary for production processes, since the color and polarizations have been summed over in the 
NRQCD matrix elements[9j. As the quark pair is in color octet, Ncolor = — 1. For different partial waves 

J = 0,1, 2, Nj = 1, 3, 5 respectively. 

As we stated before, k±_ has been set to zero in thus is independent of k±_. In eq.(j^, the integration over 
k± can be done independently. That is, 

J (fk^d^^^Pi^ + kj_)[-g>^G{z,kl) + H{z,k±)] 

= - g^j^C[wG{k^,hA_)G{zAVi\ - + 2.h‘\_hY)C[wH{k^,hl^H{zA'^_% 

C[u;(A:_L,/i_L)/(z,fci)] = J cfk±S^'^'>{k± + Pi±)w{k±,h±)f{z,k±), 


( 11 ) 
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where 



satisfies h\ = —1 and 


WG{k^,hi_) = 1, wg(fc_L,h±)= 2M2 (12) 

^ TT 

Substituting eq. 0 into eq.Q, one can see the hard part just depends on two tensors 

= -(fff + (13) 

From this structure, the two tensors are just two different gluon helicity amplitudes. Further simplification can be 
obtained by decomposing the tensor into structure functions. This is very easy to realize, because the independent 
momenta are just q, P 2 , /i_l, and only h± is transverse. Moreover, 75 in the partial wave projections always appears 
in pair, so there will be no e-tensor in and To make the expressions simpler, we put and 

together to form a two dimensional vector, i.e., , M.^}- The same rule also applies to Fi and Di 

defined later. The decomposition leads to six independent structure functions as follows: 




^2 


F^{P^q- + P^q^^) + F^{P^q'' - P^q>^), 


(14) 


where P 2 = P 2 — Z 2<1 and = e'^h±^, which satisfy 


p2-q = 0, gf = 


(15) 


From QED gauge invariance = 0, one has F^ = Fq = 0, which is confirmed by our calculation. 

Thus only four structure functions are nontrivial. It is clear that Fi’s do not depend on the scattering angles (0, ijj, 4>). 
Thus all types of angular distribution can be obtained by contracting with leptonic tensor L^''. The resulted 

cross section is 


da 

dz2dQ.dzid?- Pi_L 



- zt) 

Q4 (l-r2)2M^ 


^ CIivkIk] 

K=G,H 



1 

2 


sin^ 9D^ 


+ 


1 

2 


sin^ 9 cos 2(j)D 



(16) 


with Jg = G{z,k'^), fn = F[{z,k\) and wg,h their corresponding weights. So there are only three independent 
angular distributions for this process. Especially, there is a cos 2(j) azimuthal asymmetry, which will be an important 
signal for linear fragmentation gluon. The coefficients Di are the superpositions of Fi, 

Di = F3 - 1 - F4 — 2 Fi, D2 = F4 — F3, II3 = F3 -I- F4 -I- 2F2. ( 17 ) 


NRQCD matrix elements are contained in these coefficients. To see the contribution of different partial waves, we 
write out the matrix elements explicitly. 


A = A('^o)(08(^5o)) + A('5i)(08("A)) + ^D,CPG{OsCPi)) + ^ D,i^Pj)^{Os{^Pj)), i = 1,2,3. (18) 

j=o 


IV. RESULT 

As we can see, the analysis in last section is independent of the details of Eeynman diagrams. Since in TMD 
factorization the hard coefficients just receive contributions from the virtual correction [5], the formalism also applies 
to higher order corrections in as- At tree level the calculation is very simple. There are only two Eeynman diagrams, 
as shown in Figj^ 

After the calculation, we find that the partial waves ^Si and ^Pi have no contribution to all Di’s. The linearly 
polarized gluon fragmentation function only contributes to the cos 2(j) angular distribution, while unpolarized gluon 
fragmentation function does not contribute to this angular distribution at all, that is, 

F^ = 0, Ff = F3"=0. 


(19) 







5 




FIG. 3. The two Feynman diagrams for the hard subprocess at tree level. 



'5o 3p^ 3p^ 


3g2r2(T2 - l)/8 (1-3 t2)2/2 3 &P + 1 

3gV2(T3-l)/8 (1-3 t 3)2/2 -3(2r2-l) Gr'* - -M 

-3gV2(r2 - l)/8 (l-3r2)2/2 -3 1 


TABLE I. The hard coefficients for different partial waves. 


These two equations hold for all partial waves. Based on these facts, our final cross section can be written as 


da 


dz 2 d cos 9d?Pi_L 


(1 _ |c[wgG] ^ sin^ dUf ^ +C[whH]D^ Q sin^ 6» cos 2(/) 


( 20 ) 


and 


, 16 ( fV 2_;^)2 

^3M3(l-r2)2 

16(jV,^ - 1 )^ 
3M|(1-t2)2 


-1 


9 s 


J=0 

^ 2 
d!^ ^d^Cs^){OsCs^)) + Y^d^{^Pj) — {Os{^p.,)). 


J=0 


( 21 ) 


This is our main result. The corresponding coefficients Di can be found in Table|^ Notice that in eq.(20l zi and '0 
have been integrated over since zi is just contained in a delta function and the hard coefficients are "0 independent. 
In experiment this cross section is easier to measure. Now it is clear that the gluon TMD fragmentation functions 
can be extracted from the three independent angular distributions. By fitting the experiment data, the constraint on 
the involved four color octet matrix elements can also be obtained, which are not determined very well in literature) 
see [12] and references therein). Two features of our tree level result should be stressed here. One is about the Z 2 
dependence of the hard coefficient. As can be seen from Table|T| these hard coefficients are Z 2 independent. This 
feature will remain to higher order corrections, because the hard subprocess just depends on Pi and q, which are 
Z 2 independent. Another feature is the threshold enhancement which appears when —>■ 1. This can be easily 
understood. When the fragmenting gluon becomes soft, the intermediate heavy quark propagator, as shown in Figj^ 
is approaching to the mass shell, and this causes a factor l/Pf. Similar enhancement also appears in phase space 
integration, i.e., dPf oc Q^dzi/Pf. Since near threshold Pf is a small quantity proportional to Vl — r2, this factor 
results in an enhancement. This feature will also remain to higher order corrections. Since our knowledge about gluon 
TMDFFs and the four involved NRQCD matrix elements is very limited, here we cease to give a numerical estimate 
to the cross section. 

Before extracting these nonperturbative quantities from experiments, the most urgent task is to give a clear ex¬ 
amination for the TMD factorization for this process at least to one-loop level, since the kinematics is so simple. At 
one-loop level, two types of divergence will appear. One is caused by the collinear gluon connected to the fragmenting 
gluon, the other is caused by the soft gluon connected to the heavy quark. For the former, collinear power counting 
works, the divergence will be absorbed into the fragmentation functions by using Ward identities since all external 
particles of the hard subprocess are on-shell. For the latter, the situation is a little more complicated. There are two 
cases. First, an octet heavy quark pair from the hard interaction may emit a real soft gluon to transmit to a color 
singlet quark pair. Up to the power of v we considered here, the singlet can only be QQi^Si), i.e. an S-wave singlet. 
In this case the soft divergence is cancelled out after one sums up the divergences from heavy quark and antiquarkj?]. 
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Thus color singlet will not affect the factorization even at one-loop level. Second, after emitting a soft gluon, the color 
octet heavy quark pair may be still in a color octet. In this case, the soft divergences cannot be cancelled out by 
summing up all relevant diagrams, since we have detected the momenta of final particles. So an additional soft factor 
is required. Besides the soft divergences, to higher order of as the qq channel, i.e., 7 * qq ^ qq + QQ may break 
TMD factorization. The heavy quark pair is generated by a fragmenting gluon, and the final state pion is generated 
by a fragmenting quark. In this case, the fragmenting gluon has a off-shellness of Mj ~ '^Mq, then the contribution 
from this channel will be suppressed by Pi±/Mj, and should be ignored since both Mj and are hard scales of 
the same order in our interested kinematical region. On the other hand, the gluon TMD fragmentation functions 
may also appear in some SIDIS processes, such as e + P^A + B + X |13IfT5] where A,B are two hadrons almost 
back-to-back. It is interesting to see whether the TMD fragmentation functions there are the same as these we used 
here. The detailed discussion of these issues is beyond the scope of this paper and will be put into a future paper. 


V. SUMMARY 

In this paper we analyzed the back-to-back J/ip tv associated production at e“'"e“ colliders, making use of TMD 
factorization and NRQCD. The hadron frame where the final pion is along -|-z-axis is convenient for the analysis. 
In this frame, we find three independent angular distributions which are sensitive to the unpolarized and linearly 
polarized gluon TMD fragmentation functions. Especially, the linearly polarized gluon fragmentation function will 
contribute to the azimuthal asymmetry. At tree level, NRQCD matrix elements (08(^S'o)) and {Osi^Pj)) contribute 
to these angular distributions. Their information can also be extracted by fitting the data at e"'"e“ colliders. 
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